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METHODS OF PREDICTING HELICOPTER STABILITY 
By Robert J. Tapscott and F. B. Gustafson 


SUMMARY 


Some of the methods of predicting rotor stability derivatives have 
been reviewed. The methods by which these rotor derivatives are employed 
to estimate helicopter stability characteristics have been summarized. 

It is concluded that, although these methods are not always feasible 
for predicting absolute values of the stability of the helicopter, the 
effects on stability of changes in individual derivatives can generally 
be estimated satisfactorily. 


INTRODUCTION 


In order to predict helicopter stability - for example, to esti- 
mate theoretically whether a prospective helicopter will meet the 
flying- qualities requirements - both the applicable equations of motion 
and the necessary stability derivatives must be determined. 

The processes for handling equations of motion have been well 
established in conjunction with airplanes, and the modification of 
these procedures for helicopter use has been found a secondary problem 
in comparison with the provision of values of stability derivatives. 

The prediction of stability derivatives, in general, requires a 
knowledge of the contributions of both the rotor and the fuselage. 
Fuselage characteristics are not open to as specific an analysis as 
rotor characteristics , and preliminary estimates can be handled on the 
basis of data from previous designs and from wind-tunnel model tests. 
The fuselage seems to be subject to greater modification up to the time 
production starts; consequently, this paper is confined primarily to 
the rotor contributions to the derivatives. In the past few years, 
considerable information has been published which permits the most 
pertinent rotor derivatives to be predicted with fairly good accuracy 
where no stall is present. A general picture 
herein. 


-stV 
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SYMBOLS 
• *** 

gross weight of helicopter, 
mass of helicopter, slugs 
blade radius, ft 
blade- section chord, ft 

equivalent blade chord (on 


rotor solidity, bc e /jcR 

instantaneous blade-section pitch angle; angle between line 
of zero lift of blade section and plane perpendicular to 
rotor shaft, 0 - A-j_ cos sin ty, radians 

collective pitch; average value of 0 around azimuth, 
radians 

mass density of air, slugs/cu ft 

mass constant Of rotor blade, expresses ratio of air forces 
to mass forces, pcaR^Ip; also, climb angle, radians 

helicopter rolling moment of inertia about center of 

p 

gravity, slug-ft 

helicopter pitching moment of inertia about center of 
gravity, slug-ft^ 

helicopter yawing moment of inertia about center of 

p 

gravity, slug-ft 

difference in collective-pitch angle of front and rear 
rotors, positive when rear rotor is greater, radians; 
also increment of 0 
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AT 

AR,A(*),Aa 


difference in thrust of front and rear rotors, positive 

when thrust of rear rotor is greater, lb 

* 

definitions analogous to that for A9 


Ax d 


A 1> B 1 


5 

V 

v 

a 


a 

In- 


difference in angle of attack of front and rear rotors due 
to longitudinal swashplate tilt, positive when rear-rotor 
angle of attack is greater 

coefficients of -cos \|r and -sin i|r, respectively, in 
expression for 0; therefore lateral and longitudinal 
cyclic pitch, respectively, radians 

control motion, inches from trim 

true airspeed of helicopter along flight path, ft/sec 

sideslip velocity of helicopter in Y-direction, ft /sec 

rotor angle of attack; angle between flight path and 
plane perpendicular to axis of no feathering, positive 
when axis is inclined rearward, radians 

rotor angular velocity, radians/see 

tip-speed ratio, — — assumed equal to — 

HR m 


Up component at blade element of resultant velocity perpen- 

dicular to blade -span axis and to axis of no feathering, 
ft/sec 

Uip = Urji I OR 

blade azimuth angle (measured in direction of rotation 
from downwind position when flight-path axes are used, 
or measured in direction of rotation from position of 
blade when pointing rearward along longitudinal axis 
of helicopter when body axes are used), radians 

a slope of curve of section lift coefficient against section 

angle of attack, per radian 


L rotor lift, lb; also, rolling moment, ft-lb 

rotor thrust, component of rotor resultant force parallel 
to axis of no feathering, lb 


T 
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q 
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i 
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rotor-shaft torque, lb-ft 

T 

rotor thrust coefficient, — — 

*R 2 p(0R) 2 

rotor-shaft torque coefficient, - 

jtR^OR)^ 

blade flapping angle at particular azimuth position, 
radians; also sideslip angle, radians 

constant term in Fourier series that expresses P; there- 
fore, rotor coning angle 

coefficient of -cos \|r in expression for P; therefore, 
longitudinal tilt of rotor cone 

coefficient of -sin ^ in expression for P; therefore, 
lateral tilt of rotor cone 

projection of angle between rotor force vector and axis 
of no feathering in plane containing flight path and 
axis of no feathering 

helicopter rolling velocity, radians/sec 

helicopter pitching velocity, radians/sec 

helicopter yawing velocity, radians/sec 

angle of roll, radians 

yawing moment, lb-ft 

pitching moment, lb-ft 

height of rotor hub above helicopter center of gravity, ft 

angle between principal longitudinal axis of inertia of 
helicopter and flight path, positive nose up, radians 

Dots over symbols indicate derivatives with respect to time. 


ROTOR STABILITY DERIVATIVES 

In general, calculation of the helicopter stability derivatives 
needed for a study of helicopter characteristics depends on the knowledge 



NACA RM L5^G05 


5 


of the individual rotor derivatives shown in the following table: 


d (°T/o-) 

d(Gj/a-) 

d(Cj/a) 

d { c J/ (T ) 

a(c T /<r) 


da 

av 

aft 

d9 

aq 


dO 

aa' 

da' 

da 

aa' 

ab' 

da 

dV 

aft 

ad 

aq 

ap 

d(c Q /a) 

a(c Q /o-) 

a(c Q /cr) 

a(c Q /cr) 

d( c Q /cr) 


da 

av 

aft 

ad 

aq 


Q 

O 

o 

o 

aa 0 

O 

O 

O 

D 


da 

av 

aft 

ad 

aq 


d d| 

ao| 

aa. 

aa, 

a a. 


da 

av 

Ha 

TF 

Iq - 


ab, 

ab, 

ab, 

ab, 

ab, 


da 

"av” 

“aft" 

ad 

aq 



As shown the rotor parameters considered are Cq/cr and a', the magni- 
tude and tilt of the thrust vector, respectively; rotor torque Cq/ct; 
blade coning angle a 0 ; and bq and aq, the lateral and longitudinal 
flapping, respectively. These rotor parameters are functions of five 
independent variables: rotor angle of attack a, forward speed V, 

rotational speed ft, collective pitch 9, and pitching velocity q. The 
values of the derivatives are determined by the variations from trim of 
the rotor parameters with changes in each of the five variables. The 
total number of derivatives shown here is large; however, it is generally 
possible to reduce considerably the number needed for most applications. 
For instance, if the rotor under consideration has no flapping hinge 
offset, the flapping coefficients are not significant; thus, the number 
of derivatives to be considered is immediately reduced by about one-half. 
If necessary, however, these derivatives can be determined from the equa- 
tions of references 1 and 2. 

Of the remaining derivatives, the lift due to pitching and torque 
due to pitching generally can be neglected. For most cases, then, the 
derivatives needed have been reduced to those enclosed within the lines 
in the preceding table. All the Cq/cr, Cq/ cr, and a 1 derivatives can 
be obtained from figures or equations of references 2 and 3, as will be 
shown in the following discus si 

C 
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The thrust due to angle of attack and collective pitch are presented 
as functions of tip- speed ratio as shown in figure 1. If the solidity 
and tip-speed ratio for a given case are known, these quantities can he 
read directly.' The equations on which this figure is based and the 
processes by which they were derived are discussed in reference 3» 


The change in thrust coefficient with tip- speed ratio is presented 
in reference 3 in the form of an equation from which, when the flight 
condition is known, the rate of change of thrust coefficient with for- 
ward speed and rotational speed can be computed. 


The change in inclination of the rotor force vector due to steady 
pitching or rolling velocity has been derived in reference 2. This 
derivative is shown in figure 2 as a function of the parameter 

As the effect of tip- speed ratio p is small, the equation shown in 
figure 2 is fairly accurate below a tip- speed ratio of 0.5 for both roll 
and pitch. The rotor damping moment is determined simply by multiplying 
the quantity obtained from the figure by 18/ 7ft and by the rotor thrust 
and rotor height above the helicopter center of gravity; that is, • 


3M = da ’/dq x Tfa 18 
Sq 18/7*1 7& 


For the rotor- vector angle and rotor torque derivatives, charts such 
as those shown in figure 3 have been derived and are published in refer- 
ence 3 • These charts are given for a range of collective-pitch angles 
from 0 to 14° at 2° increments; the sample shown herein is for a collec- 
tive pitch of 8°. In these charts, the longitudinal rotor- vector 
angle a' is plotted against thrust- coefficient — solidity ratio for 
specified values of tip- speed ratio. Lines of constant power param- 
eter Cq/pCij are cross-plotted on the charts. Combinations of these 
parameters which result in angles of attack of 12° and l6° on the 
retreating blade are indicated by dashed lines; these lines, in effect, 
serve as limit lines above which account must be made of stall. By 
using slopes or differences from these charts in conjunction with other 
derivatives and some simple equations, the remaining a' and torque 
coefficient derivatives can be obtained. As an example, consider the 
derivative of a' with respect to angle of attack. From this chart, 
the rate of change of a' with thrust coefficient can be obtained at 
a given tip- speed ratio and thrust coefficient by scaling off the slope 
of the tip-speed-ratio line at the desired thrust coefficient. When this 
quantity is multiplied by the change in thrust coefficient with angle of 
attack, which has already been discussed, the result is the change in a' 
with angle of attack. The other a' and the torque coefficient deriva- 
tives can be obtained by similar procedures. These procedures are dis- 
cussed in reference 3« 
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PREDICTION OF HELICOPTER STABILITY 


By using the rotor derivatives that have been discussed, it is 
believed that the rotor contribution to the essential helicopter deriva- 
tives can be predicted. The rotor derivatives discussed axe applicable, 
in most cases, to a study of the stability characteristics of either a 
single- or tandem-rotor configuration. The difference arises in using 
the rotor derivatives to determine the helicopter derivatives for use 
in the equations of motion and in accounting for the effects of flow 
interference for a specific configuration. For most purposes, the 
lateral and longitudinal characteristics of the helicopter can be 
studied separately. Generally, equations of motion derived on the basis 
of constant forward speed are sufficient and are applicable to both 
single- and tandem-rotor configurations. 


In a study of the longitudinal characteristics of the helicopter 
(refs. 4 and 5), an important criterion is that the time history of the 
normal acceleration shall be concave downward within two seconds after 
a step input to the longitudinal control; that is, the slope of the 
normal acceleration curve shall reach its maximum value and begin to 
decrease within two seconds . In order to assist in estimating theoret- 
ically whether a prospective helicopter will meet this criterion, the 
following equations of motion were devised: 


dL n 

oq 


+ ^ A3 

de 


dL A 

— Ax - 

da 


WV 

g 


7 = 0 


dM 

dq 


1 + 


dM 

de 


A0 + 


dM 

da 


Aa + 


dM 

dB n 


AB 


x - I Y q = 0 


Aa = / q dt - A7 - ABj_ 


'0 


These equations are based on flight-path axes and derived on the assump- 
tions of constant forward speed, constant rotor speed, and constant sta- 
bility derivatives. The assumption is also made that the dynamic maneu- 
ver can be represented by a series of static conditions; hence, the rotor 
parameters are always at their equilibrium values as determined by the 
instantaneous values of angle of attack, pitch angle, tip-speed ratio, 
and pitching velocity. The first and second equations shown represent, 
respectively, equilibrium normal to the flight path and equilibrium in 
pitch. The third equation simply relates the variables of the first 
two equations to permit a solution. 
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The form as shown applies specifically to the single-rotor config- 
uration. For the tandem configuration, the pitching moment due to con- 
trol motion results primarily from differential collective pitch on the 
two rotors rather than from cyclic pitch; thus, the term in the equations 
expressing pitching moment due to control motion must he modified to 
account for this difference. 

In reference 5, these equations of motion have been solved for the 
climb angle 7, which, in turn, permits an expression to he written for 
the time history of the normal acceleration. Values of combinations of 
derivatives have been determined which, when substituted into the expres- 
sion for the time history of normal acceleration, will result in a time 
history that is concave downward at 2 seconds. These values have been 
presented in the form of a chart, as shown in figure 4, The curve in 
figure 4, therefore, indicates a boundary line separating combinations 
of significant longitudinal stability derivatives which result in satis- 
factory characteristics or unsatisfactory characteristics according to 
the criterion previously mentioned. Shown along with the theoretical 
curve are data points corresponding to five helicopter configurations. 

The derivatives for the configurations corresponding to the points shown 
were measured in flight and the adjacent number is the approximate time 
for the corresponding normal-acceleration time history to become concave 
downward. The theoretical curve is indicated to be qualitatively correct 
for both single- and tandem-rotor helicopters for separating configura- 
tions which have satisfactory maneuver stability from those which have 
unsatisfactory maneuver stability according to the criterion. This 
indication can be seen by noting that the points for which the normal- 
acceleration time history becomes concave downward in less than 2 seconds 
fall in the satisfactory region whereas, for times of more than 2 seconds, 
the points fall in the unsatisfactory region. 


Figure 4 shows an increase in the angle- of- attack stability or 
damping in pitch (that is, more negative or Mg), or increases in 

the lift-curve slope I^ to be stabilizing. Also, it can be deter- 
mined from this plot that an increase in pitching moment of inertia Iy 
can be destabilizing. Since the change in lift slope can be small only 
and is not expected to change sign, the primary changes in stability 
must be brought about by changes in the damping- in-pitch or angle-of- 
attack stability. 

In the discussion of longitudinal stability so far, the forward 
speed has been assumed constant. For the tandem configuration, however, 
the downwash effects of the front rotor acting on the rear rotor cause 
an unstable variation in pitching moment with speed, A study of the 
speed stability of the tandem was made in reference 6. In that study, 
on the basis of the available derivatives discussed previously, an 
expression for the change of stick position with speed was derived. The 
equation, along with a plot of constants based on the rotor derivatives, 
is shown in figure 5. The constants for use in the equation are presented 
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for several values of rotor solidities over a range of tip-speed ratios. 
After a flight condition is selected and the tip-speed ratio and thrust 
coefficient are thereby established, the slope of the stick-position 
variation with speed can be determined. This equation takes into account 
the effects on speed stability of center- of- gravity position, differential 
rotor speed, differential rotor radius and solidity, and "longitudinal 
dihedral. Effects of these parameters can be studied either together 
or separately. 

This analysis of tandem speed stability is of value especially in 
the study of the relative effects of changes in the various parameters 
wherein it is not necessary to know accurately the value of downwash or 
the fuselage contribution. In order to predict the absolute value of 
speed stability, the downwash must be estimated for the configuration 
and flight condition under consideration and the variation of the fuselage 
moments must be known. The results of flow-field measurements presented 
in references 7 and 8 should be useful in estimating the magnitude of 
the downwash applicable to a specific condition. 


The next important item in helicopter stability is the lateral- 
directional characteristics. If roll, yaw, and sideslip are considered 
as degrees of freedom, the equations of motion are as follows: 

IX'0 + ^ 0 + — v - ^ 6 = 0 
op ^v 06 


t «v jdN 7 . , on 
I 7 ifr + tt— ir + — v - 

Z * Sr V dv & 


Sn 




6=0 



T0 


Tty tan 7 + mv - 


cy 

cV 



For convenience in determining the derivatives, these equations are 
based on principal axes of inertia rather than on the relative wind. In 
general, the rotor contribution to the derivatives needed in these equa- 
tions can also be determined from rotor theory. Fuselage contributions 
can be predicted from wind-tunnel data and previous experience or from 
flight data where available. 


These equations have been found useful particularly in estimating 
the effects on helicopter characteristics of changes in the individual 
stability derivatives. As an example, figure 6 shows the theoretically 
predicted time histories of rolling velocity and sideslip angle of a 
tandem- rotor helicopter. The curves on the left-hand side of figure 6 
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represent the predicted time histories for the original helicopter. On 
the right-hand side of figure 6 are time histories of the same quantities 
when the effective-dihedral derivative is reduced 50 percent. From this 
figure, it is indicated that a 50-percent reduction in the effective 
dihedral would substantially improve the oscillatory characteristics of 
the helicopter. In an attempt to improve the characteristics, this means 
was tried experimentally. The results are shown in figure 7 wherein the 
experimentally measured time histories of rolling velocity and sideslip 
angle before and after the derivative change are compared. The compar- 
ison shows that, as predicted by the theory, a reduction in the effective 
dihedral improved the oscillation. Thus, in this case, the theory was 
employed successfully to indicate the course to be followed in making an 
improvement. 


CONCLUDING REMARKS 


In general, on the basis of the studies discussed herein, it will 
frequently not be feasible to predict accurately absolute magnitudes 
for the stability of the complete helicopter particularly because of the 
difficulty of predetermining final, full-scale fuselage characteristics. 
However, a first approximation can be made, and by making some compara- 
tively straightforward flight measurements of stability derivatives and 
re-employing the theory to show what modifications are needed, it appears 
feasible to handle at least those problems with which direct experience 
has been had. It appears likely that, in most cases, changes in several 
derivatives simultaneously would be necessary to achieve most efficiently 
the desirable stability characteristics, and the equations discussed 
herein should prove very useful in this type of study. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 21, 195^- • 
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VARIATION OF Cj/cr DERIVATIVES WITH RESPECT TO y. 



TIP-SPEED RATIO, y 
Figure 1 


CHART FOR DETERMINING ROTOR DAMPING 
IN PITCH AND ROLL 



Figure 2 
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CHART FOR DETERMINING DERIVATIVES OF LONGITUDINAL 
TILT AND TORQUE-COEFF. — SOLIDITY RATIO 
8 = 8 ° 



Ct /or 


Figure 5 


MANEUVER-STABILITY CHART FOR 2-SECOND CRITERION 



ANGLE-OF-ATTACK STABILITY PARAMETER 

M a Mq ( M q \2 

^ + 0 . 70 + 0 . 58 -^ + 0 . 12 ^) 

g ■ 

W V La 


Figure 4 
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TANDEM SPEED -STABILITY EQUATION 
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+ K^Aa^ + K 4 Cj 
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Figure 5 


PREDICTED EFFECT OF REDUCTION IN EFFECTIVE 
DIHEDRAL ON LATERAL OSCILLATION 




Figure 6 
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MEASURED EFFECT OF REDUCTION IN EFFECTIVE 
DIHEDRAL ON LATERAL OSCILLATION 



0 4 8 12 16 0 4 8 12 

TIME, SEC 


Figure 7 
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